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ABSTRACT 

3C 216 has a weak accretion flow luminosity, well below the Seyfertl/QSO 
dividing line, weak broad emission lines (BELs) and powerful radio lobes. As 
a consequence of the extreme properties of 3C 216, it is the most convincing 
example known of an FR II radio source that is kinetically dominated: the jet 
kinetic luminosity, Q, is larger than the total thermal luminosity (IR to X-ray) of 
the accretion flow, L^oi- Using three independent estimators for the central black 
hole mass, we find that the jet in 3C 216 is very super-Eddington, 3.3L Edd < Q < 
lOLEdd, where Q is the long term time averaged Q(t), calculated at 151 MHz. 
It is argued that 3C 216 satisfies the contemporaneous kinetically dominated 
condition, R(t) = Q(t)/ Lb i(t) > 1, either presently or in the past based on 
the rarity of L^oi > LEdd quasars. The existence of R(t) > 1 AGN is a strong 
constraint on the theory of the central engine of FRII radio sources. 

Subject headings: quasars: general - - individual (3C 216) — galaxies: jets — 
galaxies: active — accretion disks — black holes 

1. Introduction 

This article describes the extreme properties of the radio source 3C 216 and the implica- 
tions to the theory of the central engines of FR II radio sources. It has powerful lobe emission 
that implies an enormous time averaged jet kinetic luminosity (Q ~ 1.5 x 10 46 ergs/s), nar- 
row emission lines (~ 1800 km/s), and it has a very low accretion flow thermal luminosity, 
Lb i < 10 45 ergs/s (below the Seyfert 1/QSO dividing line). If it were not for the power- 
ful jet seen in a pole-on orientation, 3C 216 would be a very ordinary Seyfert 1 galaxy: 
Lboi ^ 10 44 ergs/s and L bo i/L Edd « 0.05. This object is unusual in that Seyfert 1 galaxies are 
typically radio quiet and Q in 3C 216 is an order of magnitude larger than that of any other 
jet known to emanate from a Seyfert 1 nucleus (Punsly et al 2005). The Mid-IR sample of 
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Ogle et al (2006) made it clear that many FR II narrow line radio galaxies (NLRGs) had 
no hidden quasar, but it was unclear if there was a faint central broad line AGN and how 
many (if any) of the objects had unobstructed lines of sight to the accretion disk. The blazar 
core in 3C 216 makes the orientation unambiguous, there is a direct line of sight to the BEL 
region (BLR, hereafter) and accretion disk. The Lja line is broad with FWHM = 2600 
km/sec, in an HST spectrum, Wills et al. (1995), and anything over 1500 km/sec is "out 
of family" with even the most luminous, high redshift NLRGs (Patrick McCarthy private 
communication 2006). Thus, 3C 216 is a BEL AGN and is clearly distinct from a typical 
obscured FRII NLRG (Antonucci 1993). These properties allow one to assert that the stan- 
dard astronomical estimators indicate a very small central black hole mass M bh ~ 1O 7 M 
and the jet has a time averaged Eddington luminosity, 3.3 < Q Edd = Q / L Edd < 10. 

It is not known how powerful an FRII jet can be relative to the thermal radiative 
luminosity from accretion, L^oi- Understanding the limits of jet power can help reveal the 
physical nature of the quasar central engine. The primary conundrum posed by such a task 
is that Lb i and Q need to be estimated contemporaneously and there is no reliable method 
for measuring the instantaneous Q(t) (see section 2). The only reliable estimates that we 
have are for Q. Some R = Q/Lm > 1 sources were found in Punsly and Tingay (2006). 
Unfortunately the Q estimate is not contemporaneous with the L^i data, so one can not say 
if these sources presently satisfy or ever satisfied R(t) = Q(t) / L bo i(t) > 1, (the energy flux 
emerging from the AGN in the jet exceeds that of the accretion disk radiation at the time, 
t). The fact that 3C 216 is so extreme, Q > 3.3LEdd, allows us to overcome the limitations 
of using Q and deduce that episodes of R(t) > 1 exist in FR II radio sources (in this paper, 
we adopt the following cosmological parameters: H =70 km/s/Mpc, = 0.7 and Vt m = 0.3 
in the calculations). There is no evidence that any quasar accretion flow can create an equal 
level of Lboi for the lifetime of a radio source, or even a short period of time for that matter. 
Thus, R(t) > 1, either presently or in the past based on the rarity of Lj, i > Lem quasars. 
This makes the most convincing argument to date for a kinetically dominated jet, R(t) > 1, 
in an FR II radio source. 



2. Estimating The Time Averaged Kinetic Luminosity 

Even though one can use the jet emission from the parsec scale radio core to estimate, 
Q more contemporaneously with the accretion flow emission as in Celotti et al (1997), such 
estimates are prone to be very inaccurate. A review of these issues can be found in Punsly 
and Tingay (2005) and a published example in which the Q is apparently miscalculated using 
radio core properties by three orders of magnitude is discussed explicitly. Thus, the most 
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accurate estimates of Q should rely on isotropic properties. An isotropic method that allows 
one to convert 151 MHz flux densities, F 151 (measured in Jy), into estimates of Q (measured 
in ergs/s), was developed in Willott et al. (1999); Blundell and Rawlings (2000), the result 
is captured by the formula derived in Punsly (2005): 

Q w 1.1 x 10 45 [(1 + z) 1+a Z 2 F 151 ] % ergs/s = 1.51 x 10 46 ergs/sec , (2-1) 
Z = 3.31 - (3.65) 

x (J(l + zf - 0.203(1 + zf + 0.749(1 + zf + 0.444(1 + z) + 0.205] ~ a125 ) , (2-2) 

where F 151 is the total optically thin flux density from the lobes (i.e., no contribution 
from Doppler boosted jet or the radio core). We define the radio spectral index, a, 
as F v oc v~ a . Alternatively, one can also use the independently derived isotropic estimator 
from Punsly (2005) 

Q w 5.7 x 10 44 (1 + z) 1+a Z 2 F 151 ergs/s = 1.41 x 10 46 ergs/s , a = 1.09 . (2-3) 

In order to utilize the LHS of these equations, one must extricate the diffuse lobe 
emission from the Doppler boosted core of 3C 216. Fortunately, there are many high dynamic 
range radio maps of 3C 216 (Fejes et al. 1992; Barthel et al. 1988; Pearson et al. 1985). The 
408 MHz MERLIN map of Fejes et al. (1992) shows a very luminous halo enveloping the 
core. The blazar-like, variable, flat spectrum radio core surrounded by a diffuse halo is the 
classic configuration of an FR II AGN viewed close to the jet axis in the "standard model" 
(Antonucci 1993). Using the numerous photometric data points in NED and the core fluxes 
from Fejes et al. (1992); Barthel et al. (1988); Pearson et al. (1985), one can fit the radio 
spectrum from 86 MHz to 37 GHz very accurately with a simple two component model. 
The core is variable and flat spectrum with an average flux density of about 0.9 Jy up to 
37 GHz, where the spectrum starts to turnover. The other component is very steep with 
a 1.09. This is the halo emission which achieves a flux density of 26.4 mJy at 151 GHz. 
The estimators (2.1) and (2.3) are most accurate for steep lobe emission as in 3C 216 (Willott 
et al. 1999; Punsly 2005). 

3. Estimating the Bolometric Luminosity 

The total bolometric luminosity of the accretion flow, L^, is the thermal emission 
from the accretion flow, including any radiation in broad emission lines from photo-ionized 
gas, from IR to X-ray. To estimate, L^, we use the composite SED in table 2 of Punsly 
and Tingay (2006), which see for details. This spectrum, in combination with the BELs, 
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represents the "typical" radiative signature of a strong accretion flow onto a black hole. This 
signature is empirical and it is independent of all theoretical models of the accretion flow. If 
£(y)obs is the observed spectral luminosity at the AGN rest frame frequency, u, then L bot is 
estimated as 

L bol = 1.35 vL( f\ ohs x 10 46 ergs/sec < 1.01 x 10 45 erg/s , (3-1) 

where L(z/) com is the spectral luminosity from the composite SED. 

The problem with applying the LHS of (3.1) to 3C 216 is that one must subtract off the 
high frequency tail of the synchrotron jet. The optical component is prodigious compared 
to the accretion disk. The optical emission in 3C 216 is marked by extreme variability, 21% 
optical polarization, Impey et al. (1991), and a steep power law in the optical, all confirming 
that the jet emission dominates. The best way to dig through this powerful jet emission down 
to the hidden accretion disk luminosity is to observe the object many times and hopefully 
measure 3C 216 in a low state. A low state was observed in Gelderman and Whittle (1994), 
a flux density of 1.7 x lCT 17 ergs/s/A at 4190 A in the QSO frame. This result inserted into 
the LHS of (3.1), yields the numerical estimate on the RHS of (3.1). However, the optical 
spectral index is very steep, a = 2.65, so it is likely that the the underlying continuum is 
still hidden by the jet and this is verified by the BEL strengths discussed below. 

Another method that is commonly used for estimating L bo i in blazars is to compare the 
line strengths to a composite SED (Celotti et al 1997; Wang et al. 2004). From Lawrence et al. 
(1996), the line strengths for Mgll and H/3 are 8.7x 10 41 erg/s and 7.6 x 10 41 erg/s, respectively. 
Since the line strengths are so weak, subtraction of the narrow line component from H/3 to get 
the broad component, H[3bc, is essential for accuracy. We employ the prescription commonly 
used, subtract 0.1 of the [OIII] 5007 narrow line (which is the typical ratio of line strengths 
in NLRGs) from the H/3 (McCarthy 1993; Marziani et al 2006). The Rp BC line strength is 
4.6 x 10 41 erg/s and the F(R(3 BC ) = 1685 km/s (where F() means " FWHM of , hereafter). 
Using the estimators in eqn(l) of Wang et al. (2004), L bo i ps 1.2 x 10 44 erg/s. There is no 
similar standard prescription for the narrow line subtraction in Mgll (McClure and Dunlop 
2004). The Mgll line, without the narrow line subtraction, can still provide a useful upper 
bound, L bo i < 1.4 x 10 44 erg/s (Wang et al. 2004). Note that the BEL estimates of L bo i and 
(2.1)- (2.3), imply that R = Q/L bol > 100! 

4. Estimating the Black Hole Mass 

We estimate M b h by two independent methods:, first we use the continuum flux den- 
sity and the FWHM of the low ionization BELs and secondly from the host galaxy bulge 
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luminosity - mass relation. The estimate, (4.1), is derived from the virial assumption in 
Vestergaard and Peterson (2006) applied to a sample of 32 reverberation measurements of 
AGN that provide the velocity (~ FWHM) and the radius of the BLR (time-delay times c). 
A best fit BLR radius - continuum luminosity curve is created for these AGN. Using this 
relation, La (5100 A) of a single epoch observation can be used as a surrogate for the radius 
of the BLR. Inserting the continuum flux from Gelderman and Whittle (1994) into eqn(5) 
of Vestergaard and Peterson (2006) yields, 

The assumption behind this method of taking the FWHM and La(5100A) from different 
observations is the following. La(5100A) in (4.1) represents the ionizing continuum flux. 
This flux originates in the accretion disk, since the optical flux from the jet is beamed 
orthogonal to and away from the low ionization BEL gas (Wills and Browne 1986). The 
accretion disk La (5100 A) is assumed to be fairly constant on the time scale of years and it 
is the masking radiation of the synchrotron jet that is changing in intensity. The low states 
gives the most accurate upper bound on La(5100A). The inequality sign in (4.1) is there 
because we don't know if the observed spectrum has some contribution from the jet that 
increases La (5 100 A) above the actual ionizing source value. This is a very low mass for an 
FRII radio source so it is fortunate that we have a second mass estimate for verification. 
Using F(Mg II) = 1790 km/s from Lawrence et al. (1996) and the continuum flux density at 
3000 A (extrapolated from the 4190A value in Gelderman and Whittle (1994), with a = 0.7 
as in the SED of Punsly and Tingay (2006)) in the estimator of Kong et al. (2006) provides 
a second estimate: 

The results in (4.1) and (4.2) are consistent indicating that the virial estimates are tracking 
the gravitational potential. 

The most obvious source of error in the virial mass estimates is the effect of near pole- 
on orientation. We quantify the orientation correction by defining a sample of extremely 
superluminal quasars with a projected apparent velocity on the plane of the sky of f3 a > 10 
(viewed virtually pole-on) and comparing their FWHM to "average" off-axis lines of sight, 
steep spectrum, radio loud quasars (SSQSRs). The f3 a > 10 quasars have mean F(R(3) and 
F(Mg II) of 3955 ± 1364 km/s, and 4125 ± 1997 km/s, respectively and the SSQSRs have 
mean FWHM of 5772 ± 3533 and 6398 ± 2668 km/s (Punsly 2007). If 3C 216 is typical of 
the p a > 10 extreme blazar population then the factor which Mbh has been underestimated 
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in (4.1) and (4.2) due to orientation effects can be approximated by 0.47 = (3955/5772) 2 
for H/3 and 0.42 = (4125/6398) 2 for Mgll. This analysis of f3 a > 10 quasars is in line with 
the orientation correction factor for flat spectrum quasars (FSQSRs) of 0.5 determined in 
Oshlack et al (2002) and the inference from the sample of Jarvis and McClure (2006) is a 
correction, (F(FSQSRs)/F(SSQSRs)) 2 = (4990/6464) 2 = 0.6. 

In the last section, it was noted that the line strengths and the continuum spectral index 
indicate that L^oi is overestimated in (3.1) by a factor ~ 8. Thus the effects of overestimating 
Lboi and the pole-on orientation almost cancel out in (4.1) and (4.2). The estimates corrected 
for orientation and and continuum luminosity are: 

M bh = (j^-yg)2-6 x 1O 7 M = 2.0 x 1O 7 M , H/3 (4-3) 
M bh = ( q^o.ss ) 1 - 7 x 10?M © = L2 x 1Q7M » M S H > ( 4 " 4 ) 

The estimates in (4.3) and (4.4) are consistent with each other, but it is best to verify 
the result with a completely independent method such as one based on the galactic bulge 
luminosity. The HST image has been given in deVries et al. (1997); Lehnert et al. (1999). 
This source is not trivial. The high resolution contour map in deVries et al. (1997) shows a 
strong optical jet directed to the southeast that is coincident with the VLBI jet, Barthel et 
al. (1988); Fejes et al. (1992), representing > 55% of the extended optical flux. There is a 
region of diffuse extended luminosity that is coincident with the northern radio lobe in the 
1.6 GHz MERLIN map of Fejes et al. (1992). To the east, there is excess emission enveloping 
the jet that is likely to be narrow line emission (McCarthy 1993). The host galaxy bulge 
represents only 12% of the total extended optical flux which equates to a k corrected value 
of Mr —21.3 (Fukujita et al. 1995). Inserting Mr = —21.3 into the best fit estimator 
derived from 72 AGN in table 3 of McClure and Dunlop (2002), yields M bh = 3.5 x 10 7 M o 
in close agreement with (4.3) and (4.4). 



5. A Super Eddington Jet 

The M b h estimates combined with the Q estimates in section 2, imply, 3.3 < Qem < 10- 
A value of Qem > 3.3 suggests that even if the jet central engine is currently in a low 
state then at some epoch during the lifetime of the source it must have had R(t) > 1. 
For example, if 3C 216 is not now nor has ever been kinetically dominated then the on 
average L bo i > 3.3L Edd for the lifetime of the source. However, since Q = (1/T) J Q T Q(t)dt, 
the peak values of the instantaneous Q(t), max t [Q(t)] > 3.3L Edd . Thus, one expects that 
max t [Lboiif)] would have to exceed 3.3Lem by a significant amount at certain epochs in 
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order for the 3C 216 to have always been in an R(t) < 1 state. This is inconsistent with 
the magnitudes of the peaks in the duty cycle of quasar L bo i(t)/ L Edd based on our current 
knowledge, as discussed below. 

Applying (4.1) to a sample of 232 broad line, F(H/3) > 2000km/s, AGN (138 radio loud 
and 94 radio quiet) that have been used in published estimates of M b h in Vestergaard and 
Peterson (2006); Gu et al. (2001); Baskin and Laor (2004); McClure and Dunlop (2001); 
Oshlack et al (2002); Wang et al. (2004); Punsly and Tingay (2005) reveal no L bol > L Edd 
sources. In regards to narrow line AGN, it was noted in McClure and Jarvis (2002) that 
many have broad UV lines which indicates that H/5 is not always an accurate tracker of 
the gravitational potential. This claim is borne out by the 13 narrow line objects with 
L bol > L Edd in the combined, broad and narrow line, sample of 268 AGN from the references 
above. Archival UV emission lines are available for 9 of the 13 objects. All 9 show broad 
UV lines. When F(MgII) is inserted into to (4.2) and F(CIV) is inserted into eqn(7) of 
Vestergaard and Peterson (2006), these sources are transformed into sub-Eddington AGN. 
The best available data indicates that L bo i > L Edd is a rare or nonexistent phase in the duty 
cycle of a quasar activity. The most likely explanation of the Q Edd ~ 5 estimates is that 3C 
216 experienced epochs in which R(t) ~ 5 — 10 as opposed to the alternative explanation 
that it had protracted phases in which L bo i(t)/L Edd ~ 5 — 10. 



6. Discussion 

The discussion of 3C 216 demonstrates that FR II sources seem to exist in a state of 
R(t) > 1. The Q Edd ~ 5 jet in 3C 216 must satisfy max t [R{t)] > 1, or else the accretion disk 
would episodically be in an un-physical state of L bo i(t)/L Edd ^> 1. This results is important 
because it reflects on the physical nature of the central engine. 

As a consequence of the near independence of the UV spectrum on the radio state, 
deVries et al. (2006); Corbin and Francis (1994), it has been argued that the black hole and 
not the accretion disk is the power source for FR II jets (Semenov et al. 2004). However, a 
large scale magnetic flux trapped in the vortex inside the accretion disk can tap the black hole 
rotational energy as seen in the MHD numerical simulations of Hawley and Krolik (2006); 
McKinney (2005); McKinney and Gammie (2004); Semenov et al. (2004). Yet, it is not a 
trivial matter to obtain the R > 1 condition. The simulations of Hawley and Krolik (2006); 
McKinney (2005); McKinney and Gammie (2004) are very similar and rely on the magnetic 
flux threading the event horizon of the black hole to source the jet. These simulations are 
hard pressed to achieve R > 0.1 for accretion disk parameters representative of quasars 
(i.e., an accretion disk radiating near its maximum thin disk efficiency, e > 0.2, Elvis et 



- 8- 



al. (2002)). This derives from the fact that the simulations of Hawley and Krolik (2006); 
McKinney (2005); McKinney and Gammie (2004) do not produce a net extraction of black 
hole energy, the black hole actually absorbs energy during the simulation. By far the largest 
R values occur in the new simulations of Hawley and Krolik (2006) that are parameterized by 
the black hole spin parameter, a/M = 0.99. For accretion disks with quasar-like parameters 
(see table 4 of Hawley and Krolik (2006)), R = 0.8. However, this is an overestimate of R 
that is driven by the numerical resolution of the model (McKinney private communication 
2006). Increasing the resolution in the numerical mesh near the equatorial plane allows 
short wavelength MRI modes to propagate in the accretion disk, which in turn increases 
the accretion rate, so R < 0.5 (McKinney and Gammie 2004; McKinney 2005). Even so, 
obtaining R > 0.1, in Hawley and Krolik (2006); McKinney (2005); McKinney and Gammie 
(2004) requires a very massive enveloping wind to collimate the jet. In the a/M = 0.99, 
R = 0.8 model of Hawley and Krolik (2006), the wind transports over lM Q /yr outward at 
0.3 c - 0.4 c for FRII quasar parameters. Although it is never stated, the thermal energy 
content must be enormous, ~ lGeV in order to balance the energy equation. If such massive 
winds existed around FRII jets they would have certainly been prominently detected in some 
frequency band even if they radiatively cool to a more reasonable temperature. Furthermore, 
in this numerical model the collimating wind transports far more energy flux than and 
Q combined, a high price for collimation. 

Alternatively, the black hole energy extraction model of Punsly (2001) and references 
therein that was numerically realized in Semenov et al. (2004) is based on large scale magnetic 
flux that threads the equatorial plane of the ergosphere, the active region outside of the black 
hole. The solutions attain a maximum R value, R max , that was calculated for rapidly rotating 
black holes, a/M ps 0.99 in eqn (10.14) of Punsly (2001). Extending this result to lower spins 

#max ~ ll^(a/M) ia5 , 0.9 < a/M < 0.99 , (6-1) 

where e is the radiative efficiency of the accretion flow, L^i = e(dM/dt) and (dM/dt) is the 
accretion rate from the disk to the black hole. The maximum efficiency from a thin disk is 
e « 0.3 (Novikov and Thorne 1973). The exponent in (6.1) was estimated by approximating 
the numerical data displayed in figure 10.13 of Punsly (2001) in the range, 0.9 < a/M < 0.99. 
The coupling of Q to e occurs in (6.1) because the ram pressure of the accretion flow regulates 
the maximum allowable magnetic pressure and flux in the ergosphere. For objects below the 
Seyfert 1/QSO dividing line, such as 3C 216, one might expect a/M « 0.95 and from 
Novikov and Thorne (1973), for a thin disk, e = 0.19: i? max ~ 10 which is probably required 
episodically in order to attain Q Edd ~ 5. 

I am indebted to Ski Antonucci for his valuable comments that help propel the devel- 
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opment of this research. I am also thankful to Wim Devries and Chris O'Dea for sharing 
their HST results. 
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